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Hydrodynamics with Triangle Anomalies 
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We consider the hydrodynamic regime of theories with quantum anomalies for global currents. 
We show that a hitherto discarded term in the conserve current is not only allowed by symmetries, 
but is in fact required by triangle anomalies and the second law of thermodynamics. This term leads 
to a number of new effects, one of which is chiral separation in a rotating fluid at nonzero chemical 
potential. The new kinetic coefficients can be expressed, in a unique fashion, through the anomalies 
coefficients and the equation of state. We briefly discuss the relevance of this new hydrodynamic 
term for physical situations, including heavy ion collisions. 
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Introduction. — Relativistic hydrodynamics is impor- 
tant for many questions in nuclear physics, astrophysics, 
and cosmology. For example, hydrodynamic models are 
used extensively for describing the evolution of the fire- 
ball created in heavy-ion collisions. The relativistic hy- 
drodynamic equations have been proposed many years 
ago jl|, such equations describe the dynamics of an 
interacting relativistic theory at large distance and time 
scales. The hydrodynamic variables are the local veloc- 
ity u^{x) (satisfying v? = —1), the local temperature 
T(x) and chemical potential(s) ^°'{x), where the index 
a numerates the conserved charges. The hydrodynamic 
equations govern the time evolution of these variables; 
they have the form of the conservation laws d^T^^'^ = 0, 
dfj.j'^^ = 0, supplemented by the constitutive equations 
which express T'^" and in terms of u^^, T, and 
These equations are the relativistic generalization of the 
Navier-Stokes equations. 

One feature of relativistic quantum field theory that 
does not have direct counterpart in nonrelativistic 
physics is the presence of triangle anomalies [1, 0] ■ For 
currents associated with global symmetries, the anoma- 
lies do not destroy current conservations, but are re- 
flected in the three-point functions of the currents. When 
the theory put is put in external background gauge fields 
coupled to the currents, some of the currents will no 
longer be conserved. 

In this paper, we show that the presence of quantum 
triangle anomalies leads to an important modification 
of the hydrodynamic equations. In other words, in a 
hot and dense medium quantum anomalies are expressed 
macroscopically. This modification should be important 
in many physical situations, including the quark gluon 
plasma where the small masses of the u and d quarks 
can be neglected. 

In the simplest case when there is one U(l) current 
with a U(I)'^ anomaly, the constitutive equation for the 
conserved current j^^ must contains an additional term 



proportional to the vorticity [li 



aTig^-' + u^und^i^^^+^u;^, (I) 

(2) 



where n is the charge density, a is the conductivity, and 
^ is the new kinetic coefficient. 

Even in a parity-invariant theory, the vorticity-induced 
current ^lo^ is allowed by symmetries if, e.g., is a chiral 
current. This term contains only one spatial derivative, 
and its effect is as important as those of viscosity or diffu- 
sion. Nevertheless, this term has been ignored so far. In 
fact, if one follows the standard textbook derivation 0], 
the new term seems to be disallowed by the existence of 
an entropy current with manifestly positive divergence, 
required by the second law of thermodynamics. 

Recently, however, calculations using the techniques of 
gauge/gravity duality @, @il3| within a particular model 
{J\f — 4 super- Yang-Mills plasma with an R-charge den- 
sity) gives a nonzero value for ^ [1, [^, . This indicates 
that the problem with the entropy current must be cir- 
cumvented in some way. 

In this paper we show that this new term is not 
only allowed, but is required by anomalies. Moreover, 
the parity-odd kinetic coefficient ^ is completely deter- 
mined by the anomaly coefficient C, defined through 
the divergence of the gauge-invariant current, = 
'^^LijFap, and the equation of state. 



2 ^^n 
3e + P 



(3) 



where e and P are the energy density and pressure. In 
the case of multiple U(l) conserved currents, the formulas 
are modified only slightly. Namely, Eq. ([3]) becomes 



3 t + P 



(4) 
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where a, b, c numerate the currents, C"''"^ is symmetric 
under permutations of indices and is determined from 
the anomahes, d^j'''' = -|C"^''=e'"'"^i^^^F^^. 

In the local rest frame of fluid, the new contribution 
to the current isj = ^ which means that there 

is a current directed along the vorticity. For example, in 
a volume of rotating quark matter, quarks with opposite 
helicities will move to opposite directions, a phenomenon 
which can be called hydrodynamic chiral separation. We 
note that the effect is present already at the level of first- 
order hydrodynamics. In this respect, this new term in 
the current is very unusual: it does not exist in nonrel- 
ativistic solutions of chiral molecules (e.g., sugar) in wa- 
ter, or in suspensions of chiral objects. In these cases, the 
chiral currents contains terms proportional to the third 
derivatives of velocity and terms quadratic in first deriva- 
tives of velocity The vorticity-induced current is 



hydrodynamic equations and e + P = Ts + fin, we find 



specific for relativistic quantum field theory with quan- 
tum anomalies. 

Entropy current in hydrodynamics with anomalies — 
For simplicity consider first a relativistic fluid with one 
conserved charge, with a U(l)"^ anomaly. To constrain 
the hydrodynamic equation, we turn on a slowly varying 
background gauge field coupled to the current j^. 
We take the strength of is of the same order as the 
temperature and the chemical potential, so Afj_ ~ 0(p°) 
and Ffii, ~ 0{p). As in first-order hydrodynamics, we 
keep terms of order 0{p) in the constitutive equations 
for T'^'^ and j'^ (or terms of order 0{p^) in the equations 
of motion). Note that is not dynamical. 

In the presence of an external background field the 
hydrodynamic equations obtain the form 



T 



-C^E-B. (8) 

In the standard treatment when the current is not 
anomalous, C = 0, this equation is interpreted as the 
the equation of entropy production. The first-derivative 
parts of the energy tensor and the current have the fol- 
lowing form 

T^"- = -r^P'^'^P-^'^idc^up+dpu^) - (C-|r/)P^''5 • u, (9) 
= -aTP'^'d^ (|;) + <JE^', (10) 

where P^^ = g^'^ -\- u^u'^ , and the entropy production 
rate is manifestly positive. However, in the presence of 
anomalies the last term in Eq. ([8]) can have either sign, 
and can overwhelm the other terms. Therefore, the hy- 
drodynamic equations have to be modified. 

The most general modification one can make is to add 
the following terms to the U(l) and entropy currents, 

= -aTP^^a, [pj + aE^' + ^uj^ + ^gB^, (11) 
fiA- = suf" - -f DBBf", (12) 

where ^, ^b, -D, and Db are functions of T and fi. Re- 
quiring d^s^ > for an arbitrary initial condition, and 
using the following identities which follow from the ideal 
hydrodynamic equations. 



F 



d^^f - CE^B, 



(5) 



where we have defined the electric and magnetic fields in 
the fluid rest frame, = F^'^u^, = le^'^'^Pu^Eo.f). 
The right hand sides of these equations take into account 
the fact that external field performs work on the system, 
and the anomaly. Note that the right hand sides are of 
order 0{p) or 0{p^) in our power counting, which are 
within the order hydrodynamic equations. 

The stress-energy tensor and the current are 



Tt"" (e + P)ut'u'' 



Pg 



(6) 
(7) 



e + P 



iof^id^P-nE^), 



(13) 



d^Bf" = -2uj ■ E + -^{-B ■ dP + nE ■ B), (14) 

one finds that the following four equations have to be 
satisfied 



0, 



(15) 



d^DB - ^Db - ^Bd,^ = 0, (16) 

2nD F 

2Db + - = 0, 

e + P T 

e + P T T 



(17) 
(18) 



where r^'' and v^^ are terms of order Oiji) which incorpo- 
rate, in particular, dissipative effects. Following Landau 
and Lifshitz, we can always require Up^T^'^ = u^v^ = 0. 
We find t''" and v^^ from the requirement of the exis- 
tence of an entropy current s^^ with non-negative deriva- 
tive, dfj,s^^ > 0. Transforming Uud^jT^^" + fid^J^^ using 



To proceed further, we change variables from fi, T to 
a new pair of variables, /2 = fi/T and P. From dP — 
sdT + ndfi, it is easy to derive 



dT 

ap 



T 



e + P 



dT 

dp, 



e + P 



(19) 



3 



Writing d^D = {dD/dP)d^P + {dD/Dfl)difl, and noting 
that diP and dip, can be arbitrary, as they can be consid- 
ered as initial condition on a time slice, Eq. ([15]) becomes 
two equations 

dD dD 2 

Using Eq. ITOl) . one finds that the most general solution 
to Eqs. dSni) is 

D = T'difi), ^=^{T^d{fl))p, (21) 

where d{p,) is, for now, an arbitrary function of one vari- 
able. Equation yields 

Db = TdBin), = ^ iTdB{p))p , (22) 



where dBifJ,) is another function of fi. From Eqs. pT]) 
and (fT8|) we get 

^-B(A^) = ^f^'(M)i '^b(m) - C'anomM = 0, (23) 

which can be integrated. We find 



dB{fl) = ^Cfl^ d{fl) = ^Cp'. 



So the new kinetic coefficients are 



2 



3e-t-P 



1 n^^ 



2e + P 



(24) 



(25) 



Extension to multiple charges. — It is easy to extend to 
the case of many charges. Here we consider only the U(l) 
charges that commute with each other. We denote the 
anomaly coefficients as C"''^, which is totally symmetric 
under permutation of indices and give the divergence of 
the gauge-invariant currents, 

d^ff" = c^^^E^ ■ B' . (26) 

The constitutive equations can is now 

ff" = n'^u^ + • • • + + CbB^\ (27) 

where 5° and are new transport coefficients. The 
entropy current is now modified to 

= sut" - + Duj^" + D%B''^' . (28) 

Repeating the calculations of the previous section, we 
find that 

D ^ \c'''"'T'^Ji''fJ'Jl\ D% = \c''^''Tp!'Jji\ (29) 



^ dp"- 



(:ah 



dp- 



(30) 



and, by using thermodynamic relations, one derives 
Eq. Q. 

Gravity calculation — The discussion above has been 
completely independent of details of the theory. We 
would like to check our formulas for the case when the 
kinetic coefficients can be calculated explicitly. In this 
paper we use a holographic model as a testing ground for 
our predictions. Namely, we look at the theory described 
by the following 5D action. 



4k 



LABCD 



AlFabFcd). (31) 



Here Latin indices A, B denote bulk 5D coordinates 
r,v,x,y, z, and Greek indices G {v,x,y,z} denote 
the boundary coordinates (v play the role of time on 
the boundary). The above action is a consistent trun- 
cation of type IIB supergravity Lagrangian on AdSs x 
background with a cosmological constant A = — 6 and 
the Chern-Simons parameter k = —1/(2-^/3) [13, In 
this case it describes A/" = 4 supersymmetric Yang-Mills 
theory at strong coupling, where the U(l) charge corre- 
sponds to one particular subgroup of SO (6) internal sym- 
metry. To keep the discussion general we will keep the k 
coefficient unfixed, and treat Eq. ([3T|) as the definition of 
our theory. 

The field equations corresponding to ([3T|) are 



Gab - GgAB + 2 [FacF^b + \gABF^j - 0, (32) 
VbF^^ + A^e^^'^^^FBci^iJi^ = 0, (33) 

where gAB is the 5D metric. Gab = Rab — \gABFi is the 
five dimensional Einstein tensor. The external gauge field 
and the current in the boundary theory are associated 
with the asymptotics of the near the boundary 

A,,(r,x)=A^(x)-^j,,(a;). (34) 

From Eq. p3p one derives the relationship between the 
anomalies coefficient C and k, 



C=: -- 



(35) 



These equations admit an AdS Reissner-Nordstrom 
(RN) black-brane solution. In Eddington-Finkelstcin co- 
ordinates, it is 

ds^ — 2dvdr — r^f{r, to, q)dv'^ + r^dx^ , (36) 



A 

A = -iT^dw, 
2r'^ 



where 



m q 
/(r,TO,g) = l-^ + ^ 



(37) 



(38) 



4 



The black brane is dual to a fluid at finite temperature 
T and chemical potential ji. The connection between the 
parameters of the metric and T and /i is 



2/x ^2t2 



(7 + If (37-1), 7: 



(7+1)' 



l + TT^, (39a) 



37r2T2 



^ V3 4 
The equation of state of this fluid is 



167rG5 



(39b) 



(40) 



In order to find the hydrod yna mic equations, we use 
the method developed in Ref . [lj| . We locally boost the 
RN metric and consider the boost velocity u^, as well as 
the mass and charge of the black hole, as slowly- varying 
function of the black-brane coordinates, and also turn 
on a background gauge field A^^. To zeroth order, the 
background we obtain is 

ds^ = -2u^dx^'dr + r^{P^u - fUf,u^)dx^'dx'' , (41) 



A = y^uadx'' + Al^dxt". 
2r^ ^ 



(42) 



By iteration we construct the corrections proportional to 
first derivatives 



(0) I (1) 

9AB = 9ab + 9ab 

Am — + 



(43) 



requiring the solution to be regular at the horizon. The 
solution is then expanded around the boundary; ^ and 
are read from the asymptotics of A^ near the boundary. 
As the result, we find 



is 



27rG5m ' 

a/3 (3i?* + to) qn 



(44) 
(45) 



where R is the radius of the horizon, R — ■|T(7+1). 
Equation (|44p is consistent with previous results of 
Refs. [i, 121], while Eq. (gS]) is a new result. Using the 
relationships ([M)) . one can show that the result coin- 
cides with the previous formulas (|25|) . computed using 
the equation of state (|40|) . 

Conclusion. — In this paper we show that the relativis- 
tic hydrodynamic equations have to be modified in order 
to take into account effects of anomalies. At nonzero 
chemical potentials, we find a new effect of vorticity- 
induced current. Moreover, the kinetic coefRcient char- 
acterizing this effect is completely fixed by the anomalies 
and the equations of state. 

It would be interesting to show that Eqs. ([3]) and 
(HI) is valid for all theories with gravity dual. Perhaps 



a proof similar to the one of the constancy of the en- 
tropy/viscosity density ratio 15, is possible. It is also 
interesting to understand how the new term emerges, for 
weakly coupled theories, within the kinetic framework. 
Finally, it is tempting to speculate that Eq. ^ can be 
derived by direct anomalies matching. At this moment 
we do not have a complete understanding of the micro- 
scopic origin of the vorticity-induced current. 

The effect found in this paper is a macroscopic mani- 
festation of one of the most subtle quantum effects in field 
theory. This modification of the hydrodynamic equation 
should affect the behavior of a dense and hot neutrino 
gas, or of the early Universe with large lepton chemi- 
cal potential. For heavy ion physics, one can draw a 
parallel with the "chiral magnetic effect," invoked to ex- 
plain fluctuations of charge asymmetry in noncentral col- 
lisions [i3,[il- 
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